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Abstract. Wilhelm et al. have recently called attention to the unidentified broad emission features near 117 nm in the solar 
spectrum. They discuss the observed properties of these features in detail but do not identify the source of this emission. We 
show that the broad autoionizing transitions of neutral sulfur are responsible for these emission features. Autoionizing lines 
of S I occur throughout the spectrum between Lyman alpha and the Lyman limit. Sulfur is a normal contributor to stellar 
spectra. We use non-LTE chromospheric model calculations with line data from the Kurucz 2004 S i line list to simulate the 
solar spectrum in the range 116 to 118 nm. We compare the results with SUMER disk-center observations from Curdt et 
al. and limb observations from Wilhelm et al. Our calculations generally agree with the SUMER observations of the broad 
autoionizing S i emission features, the narrow S i emission lines, and the continuum in this wavelength region, and agree with 
basic characteristics of the center-to-limb observations. In addition to modeling the average spectrum, we show that a change 
of ± 200 K in the temperature distribution causes the intensity to change by a factor of 4. This exceeds the observed intensity 
variations 1) with time in quiet regions at these wavelengths, and 2) with position from cell centers to bright network. These 
results do not seem compatible with current dynamical models that have temporal variations of 1000 K or more in the low 
chromosphere. 



1. Introduction 

This Letter is in response to the recent discussion by Wilhelm 
et al. ( 12005 > of emission features near 117 nm in the solar spec- 
trum that are much broader than observed emission lines in this 
wavelength range. These authors rule out a number of possi- 
ble explanations, such as groups of emission lines blended to- 
gether, but they show that the center-to-limb behavior of these 
features has a greater similarity to that of emission lines rather 
than to that of a background continuum, and the results in their 
Fig. 3 and Table 1 show more similarity with S i line emission 
than with line emission from C i, He i, or C iii. 

We identify these broad emission features as autoionizing 
transitions of S i. Calculations indicate that there are 47 lines 
of Si in the range 116 to 118 nm, of which 31 have dis- 
tinct wavelengths (i.e., some transitions share the same wave- 
length). All of these transitions are between known energy lev- 
els and have upper energy levels located above the S i ionization 
threshold. Six of the 31 are broad autoionization lines. Four of 
these six autoionization lines are at the wavelengths of the ob- 
served broad emission features. The remaining two are hidden 
by much stronger C iii emission near 1 17.57 nm. 

The autoionization lines are broad because they strongly 
interact with the ground state of S ii, while the narrow S i lines 
do not autoionize because their upper levels can only weakly 
ionize to higher or states of S ii. 
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Beyond identifying this wavelength coiTespondence, we 
calculate the solar spectrum in the 116 to 118 nm wavelength 
range, using a 1 -dimensional model of the average quiet-Sun 
chromospheric temperature distribution that is generally con- 
sistent with SUMER continuum observations, and with mil- 
limeter observations. We find that we can account for the basic 
observed features reasonably well, both at disk center and at 
the limb. 

In Sect. V we show the intensities calculated from a model 
of the low chromosphere with temperatures 200 K hotter than 
our average model and the intensities from a model with tem- 
peratures 200 K cooler, and show that this range of calculated 
intensities exceeds the range of both temporal and spatial vari- 
ations observed at these wavelengths on the quiet Sun. 

2. Line Data 

In 2004 Kurucz computed a line list, including line strengths, 
for Si using methods described by Kurucz J200 2). The 
line list and further details are given in the website 
|http://kurucz.h arvard.edu/atoms/1600 

His calculation used 2161 even levels in 61 configurations 
and 2270 odd levels in 61 odd configurations up to n = 16, re- 
sulting in 225605 electric dipole lines. Of these, 24722 lines 
in file GF1600.pos are between known energy levels and have 
good wavelengths. The rest of the lines have predicted wave- 
lengths. Once the lines were computed, the widths and asym- 
metries of the autoionizing lines were adjusted using Shore pa- 
rameters (Shore .1968J to approximately match the observed 
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autoionizing spectrum of Gibson et al. (I1986> . The calculations 
by Chen & Robicheaux J 19941 1 and by Altun J1992> served as a 
guide to understand the overlapping details and the background 
photoionization continuum. 

The adjusted autoionizing lines are in file GF1600.auto in 
the website. The radiative. Stark, and van der Waals damp- 
ing constants for the autoionizing lines have been replaced 
by: the FWHM Fshore, the asymmetry parameter Ashore, and 
the maximum cross-section Bshore- The profile is given by 

Ashoref + Bshoi-e/(e^ + 1) where e = 2(v - Vii„e)/rshoi-e- 

The data in the file GF1600.auto was substituted into 
GF1600.pos to make the file GF1600.sub which contains both 
autoionizing and non-autoionizing lines. The details for the 47 
S I lines between 116 and 118 nm, along with all other S i lines, 
can be found in that file. 

An earlier calculation by Fawcett (19S6) included only 4 
configurations instead of the 122 used by Kurucz. There are 
diff'erences between the Fawcett and Kurucz values of up to 
a factor of 2, probably caused by missing configuration inter- 
actions in the Fawcett calculations. 

3. The computed and observed spectrum 

Fig. [2 shows the observed disk-center intensity distribution 
for the average quiet Sun between 116 and 118 nm from the 
SUMER atlas of Curdt et al. ("2001) together with our calcu- 
lated intensities in the same absolute units (scale on the left). 
The lower part of this figure shows the S i photoionization cross 
section (scale on the right) which displays the six broad au- 
toionization lines. There is a close match in wavelength and 
shape between the calculated autoionization lines at 116.713 
and 117.055 nm and the observed broad emission features, 
and rough agreement at 117.256 and 117.846 nm. The lines 
at 1 17.500 and 1 17.600 nm are blended with C m emission. 

We include in this figure the calculated profiles of 14 nar- 
row S I emission lines having values of log gflarger than -3.0. 
(Those at 117.007, 117.091, and 117.188 nm are too weak to 
be apparent.) The strongest of these emission lines, at 116.198 
nm, has log gf- -0.906, compared with -0.603 for the autoion- 
ization line at 116.713 nm. The broad Si autoionization lines 
contribute to the S i continuum. The total continuum also in- 
cludes non-LTE contributions from hydrogen, carbon, silicon, 
and other atoms. In establishing the non-LTE populations of 
the S I levels we treat the autoionization lines only as part of 
the S I photoionization cross section and do not include these 
lines as explicit transitions. 

Figure 1 also shows the emission due to the strong C iii mul- 
tiplet with six component lines between 117.493 and 117.637 
nm, two Ci lines at 116.051 and 116.088 nm, and the Hei 
58.433 nm resonance line from the second-order spectrum 
which overlaps this mainly first-order spectrum at 1 16.867 nm. 
We took this Hei line from our calculated spectrum near 58.4 
nm, multiplied the wavelengths by 2 and the intensities by the 
factor 0.062 (determined from the 1st and 2nd order scales 
shown in Fig. 4 of Curdt et al.) and added the result to Fig. 1. 
The observed He i line has a Gaussian shape while our calcu- 
lated line has a central reversal. We are able to diminish or elim- 
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Fig. 1. The upper curves show how the calculated disk-center 
intensity compares with the SUMER observations in the 1 16 to 
118 nm range. The scale on the left gives the intensity, i.e., the 
spectral radiance. The lower curve is the S i photoionization 
cross section (right scale) that includes the six autoionization 
lines. 



inate this reversal by introducing flow velocities in the transi- 
tion region where the line center is formed. 

While the C m and He i lines are formed much higher in the 
atmosphere than the other emission lines in this wavelength 
range, we use the Cm multiplet as follows to determine the 
broadening of our calculated line profiles to compare with these 
observations. In order to match the observed blending of the 
C III components, caused by solar atmospheric motions as well 
as by instrumental broadening in these observations, the calcu- 
lated spectrum was convolved with a gaussian profile function 
having a FWHM of 0.015 nm. Lesser broadening applied to our 
calculated spectrum would give a separate emission peak just 
longward of the bright central peak, contrary to the observed 
partial blending of these two component lines. 

We also include the 14 lines of N i in this wavelength range 
having log gf> -3.0. The wavelength positions of the N i lines, 
along with those of S i, are indicated at the top of the figure. 
The four strongest N i lines ai-e at 116.745, 116.854, 116.389, 
and 116.433 nm, with log g/ values -0.675, -0.817, -1.038, 
and -1.249, respectively, according to the tables of Wiese et 
al. J1996> . The first of these appears on the red side of the 
116.713 nm broad emission feature. The second is obscured 
by the overlapping second-order He i line. The third appears in 
the blue wing of the S i 1 16.400 nm emission hne. The fourth 
N I line appears in our calculated spectrum just longward of 
the S 1 1 16.400 nm line but is not present in the observed spec- 
trum, despite having a strength 0.6 times that of the nearby N i 
116.389 nm line from the same multiplet. We cannot explain 
this discrepancy, but note that the stronger line is blended with 
the S I line while the weaker one is not. This could be resolved 
when we have included the effects of blending between S i and 
N I lines. 
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We note that while the nitrogen abundance is almost 10 
times that of sulfur, the N i lines in this wavelength range are 
much weaker because their lower energy level is an excited 
state, and the N i departures from LTE are greater than those of 
Si. 

The narrow emission feature at 117.195 nm is an uniden- 
tified line (W. Curdt, private communication), and is not the 
S I 117.188 nm line, which is too weak to appear. We defer 
a discussion of other lines and continua to a subsequent paper 
(Avrett, Fontenla, & Loeser, in preparation, hereafter AFL). 

4. Formation of the S i spectrum 

We solve the coupled statistical equilibrium and radiative trans- 
fer equations for a S i atomic model with 23 energy levels and 
111 line transitions. We must specify line strengths, line broad- 
ening parameters, collisional excitation and ionization rates, 
and photoionization cross sections. Few of these values are 
well known, but we can often determine through experimen- 
tation which parameters critically aff'ect the results and which 
do not. For example, the photoionization cross sections and the 
collisional ionization rates for the lowest S i levels largely con- 
trol the S I contribution to the continuum shortward of 120 nm. 
Emission line strengths are generally sensitive to collisional ex- 
citation rates. The S i lines considered here all have upper levels 
above the S ii threshold and are sensitive to the collisional cou- 
pling between these levels and the S ii continuum. 

The photoionization rates depend on integrations over a 
large wavelength range that includes a large number of lines, 
some strongly in emission. We cannot calculate the spectrum 
in the 1 16 to 118 nm region without solving the statistical equi- 
librium and radiative transfer equations for hydrogen. Si i, C i, 
and several other constituents in addition to S i. The AFL paper 
cited above will show results of non-LTE calculations that in- 
clude H, H , He i-ii, C i-iv, N i-iv, O i-vi, Ne i-viii, Na i-ii, Mg i- 
II, Al I-II, Si I-IV, S I-IV, Ca i-ii, and Fe i-ii, applied to the inter- 
pretation of the SUMER atlas of Curdt et al. between 67 and 
161 nm. Our non-LTE atmospheric modeling calculations use 
the Pandora computer program of Avrett & Loeser ( 2003 1. 

Our current working model of the low chromosphere is 
listed in Table 1. This model is similar to the average quiet 
Sun model C of Vernazza et al. ( 1981 1, updated by Fontenla 
et al. J1999> . but we have made adjustments to improve agree- 
ment with the observed distribution of brightness temperatures 
with wavelength in the millimeter range (see Loukitcheva et al. 
I2004> and with the continuum intensities in the SUMER atlas, 
including the continuum intensities shown here. This model 
will be revised further and presented in detail in the AFL paper 
based on more complete comparisons with SUMER observa- 
tions. 

The table lists, as functions of height (above Tsoonm = 1), 
the adopted values of the temperature and a broadening veloc- 
ity y. This broadening, or microturbulent, velocity is inferred 
from observed non-thermal doppler widths of lines formed at 
various heights, and is used not only for line broadening but 
also as a turbulent pressure velocity in determining the total hy- 
drogen number density A^h from hydrostatic equilibrium. The 



Table 1. Adopted low chromospheric model for the average 
quiet Sun 



Height(km) 


T(K) 


V(km s"') 




A' fcm"-') 


1750 


6785 


8.17 


4.17E+11 




1660 


6750 


7.64 


6.19E+11 


7.62E+10 


1580 


6710 


7.16 




7 48F+I0 


1500 


6670 


6.66 




s noF+io 


1420 


6615 


6.14 


1.94E-I-12 


8 83E+10 


1340 


6540 


5.58 


2 Q7F+12 


Q 63F-I-10 


1270 


6420 


5.06 


4.45E+12 


l.OlE+11 


1180 


6200 


4.32 


7.87E-I-12 


I.OIE+U 


1080 


5840 


3.37 


1.62E-I-13 


8.00E+10 


990 


5520 


2.61 


3.27E+13 


6.35E+10 


925 


5290 


2.27 


5.56E+13 


5.16E+10 


860 


5080 


2.01 


9.64E+13 


4.37E+10 


810 


4940 


1.77 


1.50E+14 


4.09E+10 


770 


4850 


1.58 


2.14E-I-14 


4.22E+10 


720 


4750 


1.35 


3.38E+14 


4.92E+10 


660 


4650 


1.09 


5.90E+14 


6.85E+10 



electron number density A^e is determined from the degree of 
ionization of the various elements in the calculation. 

5. Comparison of tlie broad emission and the 
continuum 

Consider the two wavelengths 116.316 and 116.713 nm. The 
first is a continuum wavelength relatively free of lines, and the 
second is centered on the strongest broad emission feature. 

Fig.Elshows, for the two wavelengths, as functions of both 
height h and monochromatic optical depth ry. the Planck func- 
tion B corresponding to the temperatures in Table 1, the contin- 
uum source function S , the mean intensity J, and the function 
dl/dh that gives the intensity contribution per unit height to the 
calculated emergent intensity / at disk center (ji - I) and near 
the limb (// - 0.2). The two emergent intensity values are in- 
dicated on the right in both cases. The emission feature has a 
larger intensity than the continuum because of its larger source 
funtion. 

Since S increases outwards, the formation region occurs at 
optical depths that are somewhat smaller than unity. The two 
panels show that the peak intensity contributions occur higher 
in the atmosphere than t = 1 for the disk-center intensity, or 
T = 0.2 for the limb intensity at yU = 0.2. 

In Fig. 13 we show the calculated center-to-limb and off- 
limb intensities for the two wavelengths. Both show brighten- 
ing toward the limb, resulting from the outwardly increasing 
temperature. The intensity in the emission feature has a max- 
imum just at, or inside, the limb, corresponding to the maxi- 
mum of S in the upper panel of Fig. 2. The continuum inten- 
sity is smaller on the disk, and exhibits a sharp increase just at 
the limb. Above the limb, the calculated intensity in the emis- 
sion feature extends to greater heights than does the continuum 
intensity. This is expected, since the opacity in the emission 
feature is greater than in the continuum. 

The limb observations shown by Wilhelm et al. in their Fig. 
5 show brightening toward the limb at both wavelengths, with 
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Fig. 2. For the two wavelengths 116.316 and 116.713 nm: the 
variation with h and with t i of B, 5, and /, and of dljdh for 
yL-\ and [x - 0.2. The emergent intensity values for /i = 1 and 
yL - 0.2 are indicated on the right for each wavelength. 

continuum intensities much smaller than those in the emission 
feature, and a maximum intensity in the emission feature that 
appears to occur just inside the limb. 

The observed intensity in the emission feature above the 
limb, however, extends further than in our calculation. Their 
Fig. 6 shows that the separation between the peak intensity of 
the emission feature just inside the limb and 0.1 times this peak 
intensity above the limb is about 5" of arc, or about 3600 km, 
which is much larger than the 2000 km or so in our calculation. 
This can be interpreted as due to the irregularities in height 
of the atmospheric layers, which are not accounted for by the 
assumed spherical symmetry in our calculations. 

Our calculations also show that the continuum intensity 
reaches a maximum at a very short distance above the limb, 
with a maximum value just above the intensity of the emission 
feature at that location. Their Fig. 5 shows extended patches of 
continuum emission above the limb that could be due to emis- 
sion from extended inhomogeneous structures that are not rep- 
resented in our 1 -dimensional modeling. The results in their 
Fig. 8 and Table 2 indicate that the continuum intensity reaches 
a maximumjust above the limb, as our calculations suggest, but 
that the continuum emission appears to remain above that of 
the emission feature at greater heights, contrary to our results. 
However, Wilhelm (private communication) points out that the 
SUMER data available in December 2004 did not allow the 
photospheric limb position to be determined to the accuracy 
best suited for this comparision. 

We have found that the calculated continuum intensity 
above the limb in this wavelength region is quite sensitive to the 
departures from LTE in H, S i. Si i, and C i in the higher layers 
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Fig. 3. Calculated center-to-limb and off-limb intensity varia- 
tions at the two wavelengths. 

of the atmosphere, i.e., to the effective scattering that we cal- 
culate in these layers due to the large number of emission lines 
overlapping the corresponding continua. We continue to study 
these effects, particularly the relative importance of scattering 
and absorption in the many emission lines that affect these pho- 
toionization rates. Further center-to-limb observations at vari- 
ous EUV wavelenghts would be very useful for such studies. 

6. Temperature Variations 

Finally we show how the calculated intensities are affected by 
higher and lower model temperatures. We calculate models 
based on chromospheric temperatures 200 K higher and 200 
K lower than in the model listed in Table 1 . The two calculated 
intensity distributions in the 116 to 118 nm range are shown 
in Fig. 13 along with the same average quiet Sun observations 
as in Fig. 1. The 200 K changes were introduced at all heights 
from the temperature minimum region into the transition re- 
gion. From hydrostatic equilibrium, increasing/decreasing the 
temperature in the minimum region causes a smaller/larger de- 
crease of A^H with height in the chromosphere, so that the re- 
sults in Fig. |4] reflect changes in density as well as tempera- 
ture. The ratio of intensities corresponding to the temperature 
changes between -200 K and -h200 K is about 4 at the contin- 
uum wavelength 116.316 nm, and is about 3 for the emission 
feature at 116.713 nm. 

The SUMER atlas of Curdt et al. gives the observed ratio 
of bright network to cell interior as a function of wavelength. 
This ratio appears to be about 1.7 at the center of the strongest 
autoionization feature, and about 2.0 at the continuum wave- 
length 1 16.316 nm. Wilhelm et al. show, in their Fig. 4, the ob- 
served variation with time (mainly 3-minute oscillations) of the 
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Fig. 4. Intensities calculated from a model with chromospheric 
temperatures 200 K higher and 200 K lower than in Table 1, 
compared with the SUMER observations of the average quiet 
Sun. 



temperature has no outward increase. The 3D simulations of 
Wedemeyer et al. (|^04} support this view that the observed 
chromospheric emission does not necessarily imply an outward 
increase in the average gas temperature but can be explained by 
the presence of substantial spatial and temporal temperature in- 
homogenities. 

However, as pointed out by Carlsson, Judge, & Wilhehn 
J1997t . simulations that do not have a persistent chromospheric 
temperature rise do not qualitatively reproduce the behavior of 
the chromospheric emission lines which are observed to be in 
emission at all times and all locations. A recent discussion of 
this issue is provided by Possum & Carlsson 2005 

Given the variation of intensity with temperature shown 
here, and given the moderate observed intensity fluctuations, 
both in continua and in lines, we conclude that models having 
a persistent outward temperature increase, but with moderate 
temporal and spatial variations, match observations better than 
current dynamical models that exhibit very large temperature 
fluctuations. Results for other wavelengths supporting this con- 
clusion will be given in subsequent papers. 

Acknowledgements. We thank Juan Fontenla, Adriaan Van 
Ballegooijen, and the referee, Klaus Wilhelm, for their comments. 



autoionization emission. The intensity has peak-to-peak varia- 
tions of about 1 .5 in the inter-network, and about 1 .3 in network 
regions. 

The calculated range of intensities coresponding to + 200 
K is much greater than the range of observed quiet-Sun inten- 
sity variations: 1) with time, as shown in Fig. 4 of Wilhelm 
et al., and 2) with position between inter-network and network 
regions, as shown in the SUMER atlas. 

These results show the temperature sensitivity of the calcu- 
lated intensities. We do not claim to have determined the tem- 
perature distribution to within + 200 K, since changing some of 
the important rates and cross sections, and changing the treat- 
ment of emission hnes involved in photoionization, can have 
large effects. 

7. Conclusions 

We have demonstrated that the broad emission features dis- 
cussed by Wilhelm et al. are the result of S i autoionization 
transitions. Using a 1 -dimensional, time-independent model to 
represent the chromosphere of the average quiet Sun, we have 
calculated the spectrum in the 116 to 118 nm band and have 
shown that the results roughly agree with the disk-center quiet- 
Sun observations from the SUMER atlas of Curdt et al. Also, 
our calculated center-to-limb variations are similar to the vari- 
ations observed by Wilhelm et al. We find that temperature 
variations of + 200 K lead to calculated intensity variations 
much greater than the temporal and spatial variations observed 
in quiet solar regions. 

With regard to such temperature variations, we note that 
Carlsson & Stein (1995 ) regard the low chromosphere in Table 
1 as wholly dynamic in nature, with temperatures varying by 
1000 K or more as individual shocks travel through the atmo- 
sphere, and with significant time intervals during which the 
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